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is not only important to the classic signal processing but also helps to disqualify several common 
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INTRODUCTION 

The goal of many scientific efforts is to predict future 
evolution of physical systems on the basis of their known 
past behavior. One example of the very limited success of 
such an activity is weather forecasts. Even if the history 
of all important parameters like temperature, pressure, 
humidity, wind velocity, etc. is known for many years 
back at almost every point on Earth, the reasonably ac- 
curate prediction of the coming weather conditions can be 
made only for several days. One may argue that Earth's 
atmosphere is especially tough system to consider due to 
its intrinsic instabilities: Even a tiny perturbation of air 
in one place can lead to huge changes of weather on a 
distant continent. In this work we will not be able to 
deal with such instable systems either. 

The other extreme is represented by very stable sys- 
tems such as celestial objects. Centuries-long obser- 
vations of the Moon and the Sun allowed ancient as- 
tronomers to predict accurately. Moon's phases, risings, 
settings and eclipses for coming millennia without any 
knowledge of gravitational forces or Kepler's laws. Such 
precision was possible because the observations of the 
system had been made over much longer period than 
the system's characteristic time scales: days, (sidereal) 
months and years. 

Would the same quality predictions be possible if we 
observed the Moon just for fifteen minutes, i.e. for time 
much shorter than the shortest characteristic timescale? 

The situation is even more interesting in quantum 
world where according to some formulations of the time- 
energy uncertainty principle it would be fundamentally 
impossible to accurately predict evolution of a quantum 
system that has been observed only for a very short pe- 
riod of time. 

In this work we show a practical way of achieving exact 
predictions of a future based on a very short history of 
a system. Our predictions will be limited only to the 
quantities evolution of which can be well described by 
finite Fourier series. This restriction is crucial. Still, 
there are many important quantities that fall into this 



category. 

Before introducing a method designed to perform such 
a task let us restate the problem in a more formal way. 
Suppose that a continuous quantity (a signal) c{t) is given 
only in a finite time interval t £ [0, T] and that it can be 
expressed in the following way 

K 

c(<) = ^4.e^'"^', (1) 
fe=i 

where K is an integer, dk is a real positive |jy amplitude 
and LUk is a real (characteristic) frequency. The quantity 
c(t) is, for simplicity, a complex function, but it can be 
made real by appropriate addition of terms with opposite 
sign frequencies ujk> — —uJk- 

Our goal is to find unknown amplitudes and fre- 
quencies ujk of given c(t) in the case when length T of 
the time interval is much smaller than the smallest char- 
acteristic timescale T <C l/u>max- The method outlined 
bellow does not require the prior knowledge of the num- 
ber of the Fourier components if in 

Mathematically, one can see that c{t) is an analytic 
function of time t and, as such, can be uniquely extended 
beyond the interval [0, T]. This means that, in principle, 
even for tiny T it is possible to get all amplitudes and fre- 
quencies from 1^ to any desired precision. Unfortunately 
it is difficult to solve this nonlinear problem analytically 
and numerical methods cannot handle continuous signals 
due to infinite number of data points. One way of get- 
ting around this problem is to take only finite number of 
points at the cost of loss of the uniqueness of the exten- 
sion. From now on we will assume that the signal c{t) is 
known only at -|- 1 equidistant time points t„ = nSt for 
n = 0, ...,N and In ^T. Q can be rewritten as a set of 
N equations 

K 

^4e-'--*" =c„, (2) 

k=l 

where c„ = c{nSt). This set has 2K real unknowns {K 
amplitudes and K frequencies) so the number of complex 
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equations A'' + 1 has to be equal or greater than K. This 
condition would have opened the possibility of existence 
of a unique solution if the equations were linear in uj^. and 
dk- It is not the case here. There will always be infinite 
number of solutions to l(2Jl if the set is self-consistent and 
no solutions otherwise. If one has an additional informa- 
tion about the range of frequencies Uk in the problem (for 
instance Moon trajectory on the sky should not involve 
frequencies higher than 1 /Hour) then the solution will be 
unique for small enough St. 

Similar problems, but for large T, are usually treated 
by a discrete Fourier transform method (DFT). That 
method assumes a grid of N equidistant frequencies and 
solves ^ only for dk as a linear set of equations. The 
spacing of the assumed frequencies is proportional to 1 /T 
and thus the method is useless for estimating values of c 
outside very short interval T . 

The more challenging task of solving (0) for both am- 
plitudes and frequencies is performed by, so called, har- 
monic inversion method. 



HARMONIC INVERSION 

The historical roots of the method can be traced back 
to the two centuries old work by Gaspard Riche (Baron de 
Prony) Q]. A numerical implementation of the original 
approach can be found in . Here we present yet another 
derivation of the algorithm. 

The key idea behind the harmonic inversion method is 
to replace the nonlinear problem |2Jl with an eigenvalue 
problem for some fictitious operator U . 

It is not necessary but very convenient to use basic for- 
malism of quantum mechanics to introduce mathemati- 
cal structure of the method. We will also benefit from 
this formalism when we turn our attention to quantum 
time-energy uncertainty principles. 

Suppose that an evolution of a normalized quantum 
state |<I>o] is generated by an unitary operator U{5t) and 
$„) = [/"((5i)|$o)- For every signal c in there ex- 
ists such an evolution operator U that the signal can be 
presented as the following autocorrelation function 



($o|$«). 



(3) 



It is enough to find eigenvalues Uk = exp{—iuJkSt) of the 
operator U (St) to find all characteristic frequencies. The 
matrix elements of U in the basis of states |<I>„) can be 
expressed in terms of c„ alone 



(4) 



where i, j — 0, N — 1. The negative indices of c in the 
equation above introduce no complication since c_„ = 
c*. To obtain K eigenvalues the dimension N of the 
matrix U must be equal or greater than K. This means 
that the number of complex signal points c„ required by 



the method exceeds the half of the number of unknowns. 
Vectors \^n) are not orthogonal, thus the eigenequation 
for matrix U takes the form 



U\uk) = UkS\uk), 



(5) 



where is a matrix of scalar products Sij = ($i|<I>j) 
Cj_i with i, j = 0, N ~ 1. The rank of the matrix S 
gives the number of Fourier components of the signal K 
for N >K. 

The harmonic inversion method consists of two stages. 
First, one numerically solves the generalized eigenvalue 
problem (|3J), in order to get all frequencies. Second, when 
all frequencies in (0) are known, a linear set of equations 
for the amplitudes dj- is solved. 

This brilliant algorithm has been developed and used 
by physicists and chemists for several years now, 0,0,0]. 
The only problem is that it fails when applied to short 
(small T) signals, 0. This limitation has been phrased 
in a form of Fourier-like uncertainty relation stating that 
the local density of frequencies that can be resolved by 
harmonic inversion must be smaller than the length T of 
the time span of the signal 

Here we claim that the applicability of the method is 
not limited by the length T of the time interval but rather 
by noise affecting the signal c{t). 

In short, solving ij^l requires calculating an inverse of 
the Hermitian scalar product matrix S. S has K positive 
and N — K zero eigenvalues. There are algebraic tech- 
niques to cope with singular matrices, see [7|. However, 
the smallest positive eigenvalue Xmin becomes very small 
for small T or large K. 



KN 



2{K-l) 



(6) 



where /(wfc) is a function of order of unity, and stands 
for the greatest frequency magnitude present in the sig- 
nal. The formula above is valid only for short signals, 
5ir <^ 1, and the expression in rectangular brackets is 
less than 1. 

Even small addition of noise to the signal may result 
in such a perturbation of Xmin that it will be impossible 
to distinguish it from, also perturbed, zero eigenvalues 
of the matrix S. When that happens harmonic inversion 
fails. The optimistic approach to © notices that Xmin 
increases very fast with increasing T. 

Assuming that the signal c{t) is corrupted by noise 
77(t), ri{t) e [-rimax,'nniax\, the ucw signal c„ = c„ -|- 
rjn has to be used to build matrices in Harmonic 
inversion will extract all frequencies if 



A. 



(7) 



which is a necessary condition assuring that all K positive 
eigenvalues of S can be found. Presence of noise leads to 
a set of perturbed frequencies Qk , which differ from ujk 



\ujk - WfelT < 



2KN^ 



An 



-Vn 



(8) 
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Notice that this is a "certainty" rather than uncertainty 
relation. The accuracy of frequencies found with the har- 
monic inversion can be made as high as needed by reduc- 
ing the ampHtude of noise rj^nax ■ This is the central result 
of this work. 

In numerical studies where the method failed for small 
T the role of noise was played by roundoff errors. Fig- 
ure n shows an exemplary application of the harmonic 
inversion to a very short signal. The values of param- 
eters in the example were deliberately chosen to expose 
the importance of the precision of the signal. 
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FIG. 1: As an example, the harmonic inversion method was 
applied to a signal with K = 10 frequencies drawn from the 
interval (0.5, 1.0), sampled at 7V-|-1 = 14 points with T = 0.01 
(upper plot). 85 digits precision, i.e. rjmax ~ 10"***, was 
used. Based on the points of the upper plot the signal is 
reconstructed (red line) with the help of harmonic inversion. 
Exact (black line) and reconstructed (red line) signals are 
indistinguishable (middle plot). In fact, the two lines start to 
differ by 1% for t > 1,000,000. The discrepancies, that are 
due to roundoff errors (noise) of the initial points, are visible 
in the bottom plot. In the example above Xmin = 4.07 x 10~^* 
and this justifies the used 85 digits precision. As a curiosity, 
the DFT algorithm applied to this signal would give just one 
frequency lu — 0, so that the prediction would be a horizontal 
line at \c{t)\ = 1. 



Coming back to the problem of determining the future 
positions of the Moon in the sky just from a short ob- 
servation we see that it is possible under one condition: 
The observation has to be very accurate. The major ob- 
stacle in achieving required accuracy would be refraction 
of the incoming moonbeams in the Earth's atmosphere, 
which results in a significant shift of the apparent Moon's 
position with respect to the actual position 0. 

The method described above can be applied to any 
signal of classical or quantum origin as long as it has the 
form of (jSJ. When applied to quantum systems, it gives 
important insight to the problem of validity of several 
formulations of the time-energy uncertainty relations. 



QUANTUM UNCERTAINTY RELATIONS 

Uncertainty principles play a central role in quantum 
mechanics. They impose constraints on the states al- 
lowed by the theory. For example, no quantum state 
can yield a product of momentum and position standard 
deviations smaller than fi/2 i.e. ApAx > h/2, where 
momentum p and position x are a pair of canonically 
conjugate operators [x,p] = ih. This also means that 
if we prepare a large number of quantum systems, all 
in the same state, and perform an exact measurement 
of position on half of them and an exact measurement of 
momentum on the other half then the spread of measured 
values must satisfy the uncertainty principle above. 

Similar in form is the relation often provided for time 
t and energy E 



AEAt > h/2, 



(9) 



which is interpreted in various ways in literature. Unlike 
energy, time is just a parameter in quantum mechanics 
and the analogy between (|5J) and other uncertainty rela- 
tions cannot be carried too far. One of the formulations 
of the time-energy uncertainty relation is presented in 
textbooks on quantum mechanics 0,^3 as follows: The 
measurement of energy of a quantum system performed 
over time At inevitably results in inaccuracy AE, so that 
(|^ is satisfied. For more ex amp les of absurd interpreta- 
tions of the uncertainty see 

It has been pointed out by Aharonov and Bohm 
that the above interpretation is wrong and one can mea- 
sure accurate energy of a quantum system in as short 
time interval as one pleases. With the use of harmonic 
inversion we can provide a simple argument supporting 
that claim. Moreover, the argument outlined bellow is 
much simpler and more general than the one used in |l2l | . 

Suppose that the state of an isolated quantum system 
spans over a finite number K of its energy eigenstates 
|fc), so that at any time t it assumes the form 



K 

E 

fe=i 



a^e 



-iujkt 



(10) 



where ak is a complex number and eigenenergy Ek — 

flLOk- 

We will show that, in principle, not only an expectation 
value of the energy but all energies E^ that contribute to 
the evolution of the system can be measured exactly no 
matter how short the measurement time is. 

First, it has been experimentally proved that one can 
actually measure a wave function ^(x) in position rep- 
resentation. The technique used in the measurement is 
known as quantum tomography The relevant theory and 
applications are reviewed in [iSj . 

Second, being able to measure ^{x) at different times 
implies that the time autocorrelation function itself 



4 



c{t) = ($(0)|$(t)) can be measured. Moreover, accuracy 
of the estimates of c{t) improves systematically with the 
increasing number of copies of the system on which such 
measurements are performed. 

The third and the last step is to use harmonic inversion 
described above to extract all eigenenergies and respec- 
tive amplitudes from c{t). 

In short, the accuracy of obtained eigenenergies Ek is 
only limited by the precision of the measured autocorre- 
lation function c{t) and this precision can be, in principle, 
as high as one needs. 



UNBREAKABLE RELATION 

Many formulations of the time-energy uncertainty 
principle were invented using intuition or dimensional 
analysis. And, as in the case discussed above, they are 
of limited applicability or simply wrong. There is, how- 
ever, one formulation that is rigorously derived from the 
quantum theory • 

The Heisenberg uncertainty relations are manifesta- 
tions of the following theorem: If A and B are two self- 
adjoint operators and a state |^') belongs simultaneously 
to the domains of A, B, Ab, bA, A? and i?^, then 

A^Ai?>i|([i,B])|, (11) 

where [l\Af = (i^) _ (i)2 and (...) = (*|...|*). The 
uncertainty above is an intrinsic feature of the state \^) 
and has nothing to do with measurement inaccuracies. 

This theorem applied to the position and momentum 
operators yields familiar relation for standard deviations 
of position and momentum. In the case of time and en- 
ergy, however, 1)11(1 results in > since time enters the 
Schrodinger equation as a parameter rather than an op- 
erator. On the other hand if A in Hll|l is replaced with 
a Hamiltonian II and |$} is not a stationary state, then 
using Heisenberg equation of motion for an incompatible 
operator B, d{B)/dt ~ i/h{[H,B]), one arrives at 

Ah41- > ^ (12) 
d{B) - 2 ^ ^ 

dt 

the uncertainty relation of energy and something that 
has dimension of time - the lifetime of the state I^P) with 
respect to the observable B. This uncertainty relation 
cannot be broken or circumvented, it holds as long as 
quantum mechanics is valid. 

There is also a time-energy relation introduced for the 
case where only one copy of a quantum system is used 



SUMMARY 

In this work we show a practical way of extracting ac- 
curate frequencies and amplitudes from a signal that is 
available only over very short period of time. The price is 
that the signal itself must be known to a very high pre- 
cision. The required precision is not achievable in real 
world experiments where observed signals are both short 
and involve many frequencies. The situation is somewhat 
better in numerical simulations where one has more con- 
trol over generated data. The role of this letter is to 
identify the source of the difficulty with a spectral de- 
composition of short signals and present a tool to perform 
such a decomposition when possible. 

The method equips us also with a powerful argument 
against some interpretations of time-energy uncertainty 
relations in quantum mechanics. 

Harmonic inversion can also be applied to signals with 
continuous spectra. In this case it provides lowest mo- 
ments of the relevant frequency distribution. 

Harmonic inversion is superior to DFT even when long 
time signals are considered. Readers interested in testing 
the method on long signals and where complex frequen- 
cies might be involved should read 

ACKNOWLEDGMENTS 

I am grateful to Jacek Dziarmaga, Krzysztof Sacha, 
Jakub Zakrzewski and George Zweig for many stimulat- 
ing discussions. Work supported by KBN grant 5 P03B 
088 21. 



[1] Negative amplitudes would only require replacing an or- 
dinary scalar product with a symmetric scalar product 
in this Letter, see 0] for details. 

[2] Baron de Prony, J. E. Polytech., 1 24 (1795). 

[3] S. L. Marple, Jr., Digital Spectral Analysis, Prentice- 
Hall, Inc. (1987). 

[4] M. R. Wall and D. Neuhauser, J. Chem. Phys. 102, 8011 
(1995) 

[5] V. A. Mandelshtam and H. S. Taylor, J. Chem. Phys. 

106, 5085 (1997). 

[6] V. A. Mandelshtam and H. S. Taylor, J. Chem. Phys. 

107, 6756 (1997). 

[7] G. H. Golub and C. F. van Loan, "Matrix computations". 
The John Hopkins University Press, Baltimore 1996. 

[8] U. Fantz, Contemp. Phys. 45 93 (2004). 

[9] L. I. Schiff, Quantum Mechanics, McGraw-Hill Compa- 
nies, 3rd ed. (1968). 
[10] A. Messiah, Quantum Mechanics, Dover Publications 
(2000). 

[11] A. Peres, "Quantum Theory: Concepts and Methods", 

(Kluwer Academic Publishers, 1995). 
[12] Y. Aharonov and D. Bohm, Phys. Rev. 122, 1649 (1961). 
[13] M. G. Raymer, Contemp. Phys. 38 343 (1997). 



5 

[14] L. Madelshtam and I. Tamm, J. Phys. (USSR) 9 249 052107 (2002). 

(1945). 

[15] Y. Aharonov, S. Massar, S. Popescu, Phys. Rev. A 66, 



